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I— I Abstract 

^~~^ We supplement a very recent paper of G. Markowsky concerned with the 

J^ expected exit times of Brownian motion from planar domains. Besides the use 

■^ of conformal mapping, we apply results from potential theory. We treat the 

case of a wedge-shaped region exactly, subsuming an example of Markowsky. 

A number of other results are presented, including for a half disc, n-grams, and 
\1J a variety of other regions. 
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Introduction 

We will be concerned with developing exact results for the expected exit time 
of Brownian motion from certain simply connected planar domains. The domains 
include polygons, not necessarily regular, and others. The results have applications 
to other areas of both pure and applied mathematics. This Introduction also shows 
our notation. 

Intuitively, the average exit time of Brownian motion gives a sense of the size of 
the domain, and the distance from the starting point to the boundary. One may 
then make connections between such probabilistic considerations and the theory of 
analytic functions [6| ITOj. 

For this, let S be the Schlicht class of analytic functions / on the unit disc D, 
normalized such that /(O) = and /'(O) = 1. By the Riemann mapping theorem, 
there is an onto map of D to simply connected domains Q. In a sense, the extremal 
functions in S are the identity function I{z) = z (the smallest) and the Koebe function 
K{z) = z/{l — zY (the largest). Earlier results [71 [10] essentially show that the 
expected exit time of Brownian motion starting at the origin of a simply connected 
domain is bounded above and below by the corresponding times for the Koebe domain 
and the unit disc, respectively. For the Koebe domain, the expected exit time is 
infinite while for the unit disc it has the value 1/2. 

Let R be an open, connected subset of M" with n > 2, and X a Brownian motion 
starting at a point x in R. Then r is the first time X leaves R: t[{jS) = inf{t > : 
Xt{uj) ^ R}. A result that is convenient for two dimensional applications is provided 



by the following [1]. 

Lemma 1. Let f{z) = Xl^o'^"'^" ^^ conforinal on D. Then for Brownian motion 

starting at /(O), and exit time r, 

Eno)[r{f{D))] = -J2K\'. 

n=l 

This result can be proved by using the optional stopping theorem for martingales. It 
can also be developed by using the Green function for the two dimensional Laplacian. 
Another approach to finding expected exit times is through potential theory. Then 
the exit time is given by the unique solution of the Poisson equation V^m = —2 in 
R with u vanishing on the boundary of the domain (Dirichlet boundary condition) 
(e.g., [TlJ). Equivalently, if R has a suitably normalized Green function G{x,y), 
the expectation of r as a function of the starting point x is given by E^[t{R)] = 
j^ G{x, y)dy. For some highly symmetric polygonal domains, including the equilateral 
triangle and the rectangle, there are known eigenfunctions and eigenvalues from which 
the Green function may be constructed. As far as the equilateral triangle, there is 
research starting with Lame [16] and ongoing [1], [TH EU [251 ES]. Other recent work 
also suggests how to develop results for other symmetric domains including a regular 
rhombus and hexagon [191 120] • For the equilateral triangle, the eigenvalues of the 
Laplacian with zero Dirichlet boundary condition are given by \mn = ^^^ ijn^ + n^ — 
mn), m,n = 0, ±1, . . . subject to: m + n is a multiple oi 3, m ^ 2n, and n ^ 2m. In 
particular, Pinsky f2M 125] investigated the multiplicity of the eigenvalues, observing 
that any eigenvalue can be written as the norm of an integer in the quadratic field 
(5("\/— 3), and the symmetry of the eigenfunctions. 



Very recently Markowsky [18] investigated expected exit times from simply con- 
nected planar domains. In particular, an example concerns the expected exit time 
from a regular polygon Um of m sides (m-gon). In this circumstance, the well known 
mapping from the unit disc to the polygon is provided by a case of the Schwarz- 
Christoffel transformation with all interior angles equal, with value vr^™"^)/'". 

The expected exit time from Um is expressed in terms of a generalized hypergeo- 
metric function pFq, 

2 
Tfl 

E^[T{Um)\ = TTfWTT-, ^ ^^^ AF^{l/m, 1/m, 2/m, 2/m; 1 + 1/m, 1 + 1/m, 1; 1), 

(1.1) 
where B{x, y) = r{x)r{y)/r{x+y) is the Beta function and F is the Gamma function. 
The generalized hypergeometric function is defined by 

pi^,^ai,...,ap,Oi,...,o,,xj Z^ (5^)^...(5j^ A;! ' '' ^ ' 

where 

, , Tia + k) 
T[a) 

is Pochhammer's symbol. This function has a branch point at x = 1. Under suit- 
able conditions on the parameters, convergence may be obtained on the unit circle. 
Amongst many properties, the function w{x) = pFq{{ai}] {bj};x) satisfies the linear 
ordinary differential equation 

q p 



oY[{0 + bj-i)-xY[{e + 

L i=i i=i 



w = 0, p < q + I, 



when no bj is a nonpositive integer. Here, the differential operator 9 = x-^. For the 
sake of brevity, we refer to several standard sources for information on the functions 
pFq and other special functions [I El [HI [131 [13 EZl EH] ■ 

Another special function that we will use is the two-variable Appell hypergeomet- 
ric function Fi, given by the series [3] 



oo oo 



F,(a; 6„ h,- c;x,y) = J2i: i^f^M^^-^^n. (1.3) 



This function has the integral representation 



T{a)T{c-a) Jq 
and satisfies a system of two partial differential equations of second order each (e.g., 
[13], p. 1054). Some reductions for Fi are known, including 

F,{a- P, /?'; 7; X, 1) = y^ ^\ y- 2Fi{a, /3; 7 - /?'; x), (1.5a) 

1 (7 - a)i (7 - p') 

and 

Fi(a;/3,/3';7;x,x)= 2Fi{a, /3 + /3'-r,x). (1.56) 

We first present results based upon the use of conformal mapping and Lemma 
1 for wedge-shaped regions, the half disc, symmetric star-shaped regions (n-grams), 
and a region bounded by circular arcs. We then give results for a variety of domains 
by providing the solution of the Poisson equation with constant source. 

Results from conformal mapping 

The following generalizes Example 5 of [15] . 
Proposition 1. Let EI = {Re z > 0} and let IF = {|arg(2;)| < TTp/2} for p < I, 



where a.Tg{z) takes values in (— 7r,7r]. Then 

oo 

EMm] = 7,Y. —-mi V 2^i'(-'^' -P; 1 - m -p; -1). (2.1) 

m=l ^ ' 

Here, EI is the infinite wedge of width vrp centered on the positive real axis, 
and H^ = H. This Proposition confirms that [3T] £^i[r(EI^)] < oo if p < 1/2, and 
^i[r(eP)] =oo ifp> 1/2. 

Proof. The conformal map from D to IP is g{z) = (1 + ^)^/(1 — 2;)^. (As illustrated 
in Figure 1.) We develop the Maclaurin expansion of this function with a double 
application of the binomial theorem. In fact, we do this slightly more generally with 
a function gg-. 

*<^' = frS = '1 + --)'£ (T) '-')"--" = (1 + --)'E ^-" 

^ ' n=0 ^ ^ n=0 



j=0 n=0 ^-^^ 

oo oo / \ / \ 

j=0 m=j 

oo m 



j=0 m=j V"' / ^ -^ ^ 



■z'" 



m=0 j=o ^-^ / ^ -^ ^ 



oo 



^a^^-, (2.2) 



m=0 



where 



<^m=y[ . 7 ^ = p. . .^-p. ^ 2Fi{-m, -q;l-m-p; -1). 2.3 



Then by applying Lemma 1 with q = p we obtain the Proposition. For any p or q, 
the first term of the 2-Ff factor in (2.1) is 1. Then p = 1/2, where the factor r(l — 2p) 
diverges, provides the critical case. [For further details on this point, see the special 
case discussed next.] 

We may easily verify that at g = a case of Markowsky is recovered. We have 

1 1^1 T^{m + p + l) 



1 y. 1 1 _ 1 y- 



2 ^ m'^ B'^{m,p) T'^{p) "^ {rn + 1)^ r^(m + 1) 
2 ^ (m + l)2 (m!)2 

m=0 V ' ^ V ^ 



^[ 2F,{p,p- 1-1)- 1] = ^- 



r(i - 2p) ^ 



(2.4) 



[T^{l-p) 

being equivalent to (3.13) in [18]. In relation to Proposition 1, Markowsky [18] wrote 
that "the Taylor series expansion for g appears to be unwieldy for arbitrary p" . 

We next obtain a hypergeometric identity. 
Proposition 2. For < p < 1/2, 

oo ^ 1 

sec(7rp) - 1 = V ^^^7 ^ 2Fi{-m, -p;l - m - p;-l). (2.5) 

Proof. We solve the Poisson equation for a wedge-shaped region of half angle a. 
Using polar coordinates, V^f = —2 reads 

1 d f dv\ 1 d'^v 
r dr \ dr J r^ dO"^ 

We make the solution ansatz v{r,6) = r^g{6), with n and g{6) to be determined. 
Since then g"{6) + n'^g{9) = —2r^^'^, we take n = 2. With a particular solution 
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Figure 1: Conformal map from unit disc onto wedge. 



gp{6) = —1/2, the solution for v has the form 



v(r, 9) =r^ [a cos 29 + b sin 29 



(2.7) 

The boundary condition v{r, 9 = ±a) = gives a = 1/(2 cos 2a) and 6 = 0, in which 
case 

r^ / nns 9ft \ 

(2.8) 



vir,9) = '-('-^^-l\. 



2 \ cos 2a 

The wedge being an unbounded domain, there is the question of the uniqueness of 
the solution. As the solution (2.8) grows no faster than |zp, it is indeed the unique 
solution [B]. 

For the point /(O) = 1, we evaluate v{r = 1,9 = 0), putting 2a = vrp. Using 
Proposition 1, Proposition 2 then follows. D 



Corollary 1. 



1 1 



1 3 

4' 4 



(2.9) 



Remark. The Green function for the wedge can be obtained explicitly via sep- 
aration of variables and Fourier series. We could then integrate over the wedge to 
determine the solution of the Poisson equation. We have avoided that prospect. 

The following gives the expected exit time from the half disc Dj^ = {\z\ < Ijlmz > 
0} from the point /(O) = -i(l - ^2). 
Proposition 3. 

[m/2] 



^/(0)[r(/P+))] = 2 5^ 



m=l 



1- 






V2 



e=o 



1/2 



(2.10) 



Proof. The conformal map of D+ onto the unit disc is given by [22] (Example 
7.35, p. 254) 



c 



z'^ + 2iz + l 



(2.11) 



^2 -2iz + l 
Therefore, the inverse mapping of the unit disc onto the half disc that we require is 



r^(^) 



{z-l) 



\1 + z - V2VI + z^]. 



(2.12) 



This mapping is illustrated in Figure 2. We then employ binomial expansion to 
develop the power series of (2.7), 



c-(--) = E 



m=0 






.2^ 



(2.13) 



Reordering the sums here and applying Lemma 1 gives the Proposition. 



D 




2 -2 




-1 - 



-2"- 



Figure 2: Conformal map from unit disc onto upper half disc. 



Remarks. The approximate numerical value of (2.10) is 0.191807. 
The other inverse mapping from (2.11), 



C\z) 



■[l + z + V2Vl + z^]. 



(2.14) 



(^-1) 

takes the unit disc onto the region outside of the half disc in the lower half plane, and 
is illustrated in Figure 3. 

Results are known for the Green function for the Dirichlet and Neumann problems 
for a half disc [5]. Therefore, it may be possible to complement Proposition 2 with 
an approach of integrating the Dirichlet Green function [5J (pp. 257-259) over the 
half disc. We note that this method gives in principle a way to find the expected exit 
time for any point within the half disc. 
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Figure 3: Conformal map from unit disc onto lower half plane exterior to half disc. 



In fact, with the reflection symmetry, we may anticipate that a difference of Green 
functions for the circle will suffice. So for the special homogeneous boundary condi- 
tion, we have 

(2.15) 



G(z,C) = -ln 



{z-0{l-zO 



[l-zO{z-0 
Introducing z = re*^ and ( = pe*"^, if the half disc has radius r^, we have the expression 



G{z,C) 



1 



In 



r"^ + p"^ — 2rp cos{6 — 
rg + ^ - 2rp cos(^ - 



-In 



r'^ + p'^ — 2rpcos{9 + q 
rg + ¥ - 2rp cos(^ + 



(2.16) 



This Green function is symmetric in z and (, as it should be. If the half disc has 
radius tq, then the solution for the expected exit time is given by 



PIT pro 

u{r,e) = -2 / G{r,p,e,(l))pdpd(j). 
Jo Jo 



(2.17) 



11 



Noting that, for instance, 

I ln(a^ -2a6cosx + 6^)(ix = 27™ln[max(|a|, |6|)], (2.18) 

it appears that these integrals can be evaluated explicitly, but are lengthy. The 
solution for a disc is reviewed in the Appendix. As well, (2.18) is verified there, and 
other integrals are demonstrated. The latter suffice to perform (2.17) first in terms 
of dilogarithms, which details will be given elsewhere [H]. Upon combination of all 
contributions, and simplication, the solution of the Poisson equation can be written 
in terms of elementary functions. 

Indeed, from our exact solution for the Poisson equation and Proposition 3 we 
have obtained the following identity. 
Corollary 2. 

m=l Y=[m/2]+l ^ ^ j 

Similarly, the Dirichlet Green function for a quarter disc may be written. We 
expect that it too can be integrated to provide an explicit solution of the Poisson 
equation. 

We next consider n-grams, non-convex star-shaped regions of 2n vertices with 
alternating interior angles vrcti and '7rQ;2, with a\ < \ < a^- The following result is 
useful since it implies bounds for the expected exit time. We let R denote the radius 
of the circle circumscribing the ra-gram, and Ry) the inradius, the radius of the largest 
inscribed circle. 
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Proposition 4. For an n-gram of alternating exterior angles nfii and 7r/i2, 

1 (1 - /i2 + l/nj \ n n J 

and 

r(i-/ii)r(i + i/n) ^ /- 1 ^1 \ ...Q,^ 

r(l — /ii + l/n) \ n n J 

Proof. We let 7r/ij = 7r(l — Q;j),-i = l, 2be the corresponding alternating exterior 
angles. Because the sum of all exterior angles is 27r, we have /Ui + /i2 = 2/n. The 
Schwarz-Christoffel mapping from the unit disc onto the n-gram has the form 

where A and -B are constants that may be included to dilate and rotate, and translate, 

respectively, the n-gram. We take Hodd = /ii > and Heven = —1^2 > 0, and Uj as 

2nth roots of unity, so that 11 7" en (-^ ^ ^j) — ^"' ^ ^ ^^^ YliLdi^ ~^j) ~ z"" + 1. Then 

we may write 

^W = / Trr^dz. (2.21) 

Jo (1 + ^ )^' 

With a change of variable, this integral could be evaluated in terms of Fi using the 
representation (1.4). Instead, we perform binomial expansion, integrate term-wise, 
and appeal to the series definition (1.3): 

^n(i+fe)+l 



■"W-EE.-KlKT) 



j=0 k=0 



k J [n{j + A;) + 1] 



V^>yp/ -|Nfc (/^2)i (/^l)fc (l/n)j+fc njj+k) 

^^ j! k\ (l + l/n),+fc 
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= zF, (-■ fi,,fi2; 1 + -; z^- -zA . (2.22) 

To find the radius and inradius, we evaluate at the preimage points z = ±1, using 
(1.5). 

Remarks. A pentagram example is given in [21] (pp. 196-197). 

When /ii = 1 — 1/n and /i2 = S/n — 1, the 2-^1 factors in (2.19) may be reduced 
in terms of ratios of Gamma functions. 

If we reorder series in (2.22) we may write the Maclaurin expansion for the map- 
ping w{z), 

^(') = E E(-l)^' ( ) ( ■ ) 7^ rv (2-23) 

m=0 j=0 \ J / \ J / \ 'I 

Then by Lemma 1 we obtain the following. 

Corollary 3. For an n-gram of alternating exterior angles vr/ii and '7r/i2, the expected 

exit time from the origin is 

[(fc-l)/n] 



Mr{jm\ = \Y.l 



2^^P 

fc=i 



S '-nT)!-"-. 



(2.24) 



For a univalent map / of the unit disc, there is a number W, independent of /, 
such that the inradius Rf{D) > ^l/'(0)|- The number lA is called the univalent or 
schlicht Bloch-Landau constant. Upper and lower bounds are known, 0.570884 < 
lA < 0.6563937 (e.g., [Bj), based upon using certain constructed domains. For the 
n-gram family, the factor /'(O) = 1. 

For the next exact result, we let L be the region bounded by the circles Ci, 
I2; — 1| = \/2 and C2, \z + l\ = \/2 and containing the origin. 
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Proposition 5. 



Eo[r{L)] 



1 1 



(2.25) 



TT 2 

Proof. We will use the following. 
Lemma 2. A conformal map from L to the unit disc \w\ < 1 is given by w = 

Proof. For points on the circular arc Ci, z = 1 + y^e*^, with 37r/4 < 9 < 57r/4, 
so that 



w 



l + y2e 



-6^^(72 + 6 



i9\ 



This gives 



w 



1 + y2(cos 6* + i sin / 



V^ + cos 9 + i sin ^ 



;i + \/2cos^)2 + 2sin^^ 



(v^ + 



cos I 



sin'^ 



1/2 



1. 



Similarly, for points on the circular arc C2, z 
Then 



4 + V2e*^ with -n/4 <9 < n/4. 



w 



V2e 



ie 



and again \w\ = 1. 

In addition, for f{z) = 2z/{l - z^), f'{z) = 2(1 + z^)/{l - z^f. In particular, /' 
does not vanish within L. With f{z) being one-to-one in this region, it follows that 
/ is conformal. To verify that that the interior of L is mapped into the interior of 
the unit disc, we may consider the vertical line segment z = it with |t| < 1. Then 
|w| = 2|t|/(l + t2) < 1. D 

From the Lemma, the inverse mapping from the unit disc to L, illustrated in 
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Figure 4, is given by 

z{w) = -[-l + Vl + w% (2.26) 



w 
and its Maclaurin expansion is 



oo 



i:(^;>"--i:(;fO»-'- (-^) 



Then by Lemma 1, the expected exit time is given by 

2 



Eo[riL)] = -y(^ ^^\, ] =---, (2.28) 



the sum being a special case of 






2 1 r / 1 1 

1+ 2^1 (-2, "25 l5^' 



(2.29) 

D 



Remark. The mapping (2.26) is a special case of 



\{l + zf - il- zy]i 

wJz) = ^\J { )- V^. 2.30 

^^ ^ (l + z)P + (l-z)P ^ ^ 

Here the circular arcs meet at an angle of np. (2.26) then corresponds to p = 1/2 
and an inessential pre-rotation z ^ iz oi the unit disc. 

Other results from potential theory 

The following solutions of the two-dimensional Poisson equation emphasize that 
exact results for average exit time may be found for all points of certain polygons 
and other domains. Specifically, we consider solutions of {d^ + dy)u{x,y) = —2 with 
zero Dirichlet boundary condition. 
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Figure 4: Mapping of unit disc to lens-shaped domain. 

Equilateral triangle. Let the triangle of side length a with origin at the center be 
oriented as in Figure 5. The solution of the Poisson equation is 



u[x. 



y) = ——{2\/3y + a){\/3y + 3x - a){\/3y - 3x - a). (3.1) 

18a 



Should the triangle be inscribed in a circle of unit radius, then by the law of cosines 
we have a = v3. Then ^(0, 0) = a^/18 = 1/6, in agreement with EqIt^Us)] = 1/6 [1]. 
Here the area of the triangle is A = a/So^/S, and the product form of the solu- 
tion (3.1) may be recognized as parallel to that of ^ given in terms of barycentric 
coordinates (Ai,A2,A3), 2v^y4AiA2A3. At least at the time, the authors of ^ were 
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unaware of the exact solution (3.1). However, the solution (3.1) has been well known 
in torsion theory for a long time (e.g., |,3Qj ) . 




V3y - 3x - a = 



V3y + 3x - a = 



2V3y + a = 



Figure 5: Equilateral triangle domain. 



A circular domain with circular indentation. We next consider the domain of 
Figure 6 with radii a and b with a > b. The solution of the Poisson equation is given 
by [30] (pp. 126-127) 



u{x,y) 



2 \ X'' + y'^ 



(3.2) 



or, in terms of polar coordinates with r^ = x^ + y'^ and 9 = tan ^{y/x), 



u{r, 9) 



l-2a- 



cos^ 



(3.3) 



Again we know the expected exit time of Brownian motion from all points of the 
domain. A special case is that from the point at the center of the larger circle, 
u{r- = a, ^ = 0) = |(a^ — IP'). We also obtain the expected result as 6 — )■ 0. 



y 


/" r"^ 
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b . 

.■'\ 


-V"^ 


\ 




"\ \ 


/ "" 




\ a\ 


/ 



Figure 6: Circular domain with circular cutout. 

Isosceles right triangle domain. The geometry, with two sides of the triangle of 
length a parallel to the coordinate axes, and the hypotenuse along the line y = —x, 
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is detailed in Figure 7. By combining homogeneous and particular solutions, the 
solution of the Poisson equation has the form 

'1 



uix.y) = - 



^{x^ + y^) -^P{x,y) 



(3.4) 



where ^/^ is a harmonic function. We let kn = {2n + l)7r/a. Then ilj{x,y) is given by 
[30] (pp. 133-134) 

ip{x,y) = ~xy+-{x+y) ^ T^ 77^ — , .^n . , ., --(sinh fc^ycos A;„x+sinhfc„xcos/i;„|/). 

2 TT'^ -^ (2n + l)-^smh(A;„a/2) 

(3.5) 
This function has been constructed so that its restriction to the sides of the triangle 
is (x^ + y^)/2; this may be verified with a certain cosine Fourier series. Hence (3.4) 
satisfies the zero Dirichlet boundary condition. A special case of (3.4) is for points 
along X = y, wherein 



,2 °° 



u{x, y = x) = -2x^ + ax ^ V" -- — , ^,„ . , ,, r— sinh fc„x cos Kx. (3.6) 

^ ^ {^^ + 1) sinh{kna/2j 



ia^Y^ (-1] 

n=0 

One may compare the diffusion application in [22] with a = 2. 

Elliptical domain. Let an ellipse be given by the equation x^/a^ + y'^/b'^ = 1. Its 
foci are at ±/ = ±Va^ — 6^ and its eccentricity e = a/1 — h"^ /a?. The solution of the 
Poisson equation with zero Dirichlet boundary condition is 

The maximum average exit time is evident, if (0, 0) = a^lp' /{a^ + 6^). For a = b = R 
we recover the known result for a circular disc. A conformal mapping of the interior 
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Figure 7: Right isosceles triangle domain. 



of an ellipse to the unit disc is known, and involves the elliptic function sn [2T] (pp. 
295-296), [23]. 

Rectangular domain. Omitting some details, we write the solution of the Poisson 
equation with constant source —2 on a 2a x 26 rectangle in the xy plane. 



2,24^ (-l)"cosha„|/ 



?i=0 



COS a„x. 



(3.8) 



where a„ = (n + l/2)7r/a. This solution follows from the use of separation of variables 
and the application of Fourier series in the x direction. Once again there are extensive 
physical realizations of this solution, including electrostatic and gravitational poten- 
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tials, the transverse deflection of a membrane, and steady, incompressible fluid flow. 
The maximum average exit time occurs at the center of the rectangle. 



m(0,0) = a^ 



32 



E 



■1)" 



Ti" '-^ (2n + 11 



rsech 



n+ -] -b 
2 a 



From (3.8) we may then supplement relations (3.21) and (3.22) of Markowsky's work 



Eo[t{U,)] = 4^3(1/4, 1/4, 1/2, 1/2; 5/4, 5/4, 1; 1) 



52(1/4,1/2) 



64 



TC^ 



oo oo 



EE 



\m+n 



with 



-^^^ (2m - l)(2n - l)[(2m - 1)^ + (2n - 1)^] ' 
Eo[r{U,)] = 4^3(1/4, 1/4, 1/2, 1/2; 5/4, 5/4, 1; 1) ^ 



(3.9) 



52(1/4,1/2) 



9 ^^-^ r\3 



Oit. cosh «„/ a/2' 



(3.10) 

for the expected exit time from the square. For the square inscribed within the unit 
disc, we have put a = b = l/V^. It may be noted that we thus have a form of 
the expected exit time for a rectangle of arbitrary dimensions and for an arbitrary 
starting point. 

The double series form, (3.9) above, (3.22) of [18], is based upon the double Fourier 
sine series solution in [15] (p. 71). It is then apparent that that representation and 
(3.10) are related via the partial fractions forms 

'7rx\ Ax s-^ 1 



tan 



2 J TT -^ [(2A;-1 



TTX \ 1 2x 



)2-x21 



2 7 nx vr -^ (x^ — A;^) ' 



(3.11a) 
(3.116) 
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and 



4>^ (-l)H2k + l) 
secM = sech x = - }_^ \(oi. , ^^2 , ^21 - (^-Hc) 



n^\{2k + ir + x^y 



fc=0 



Indeed, from (3.11a) and (3.11b) we have 



E (L-l)|(2„.-l)^ + .^l = 8 ;3 1-2 + "t|(-/4)(- + 1)1 + ta„l(./4)(i. + 1)11 • 

m=l ^ ' '-^ ' 1 

(3.12) 
For an infinite strip of width 2a, we have the special case of (3.8) u{x) = a^ — x^ . 
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Appendix: Green function solution of the Poisson equation for a disc and integration results 

Here we first use the Green function corresponding to the first term on the right 
side of (2.16) to solve the two dimensional Poisson equation with homogeneous Dirich- 
let boundary condition. Noting that by symmetry, the solution can not depend upon 
6, we have 



I pro r-n . ^.2 ^ p2 _ 2rp cos > 



pro i'2-K 

'0 Jo 
Using (2.18), one contribution is 

^r•r() j'2tt 
io 



n + 



2rp cos 1 



pdpd(f). 



pra pZtt pro 

/ / ln(r^ + P^ ~ 2rp cos (f))pdpd(f) = An ln[max(p, r)]p(ip 
Jo Jo Jo 



Att 



TO 



\nr pdp + pin p dp 

Jo Jr 

= 7c[2rl In ro + r^ — Tq] . 



(Al) 



(A.2) 



The other contribution is 

no /-Stt 



JO 



^.2p2 \ rro 

In ( Tq H ^ 2rp cos ) pdpdcj) = An j In 



r, 











max ro. 



rp 
ro 



pdp 



ro 



= An In ro / pdp = 2nrl In Tq. 
Jo 

By combining (A. 2) and (A. 3) in (A.l), we recover 



(A.3) 



uir) = -{rl-r^). 



{A.A) 



n 
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Lemma 3. We have the expansion 



ln(r^+p^— 2rpcos(6'=l=0)) = 2 In max(r, p) — 2 ^ ^ — j 

k=l ^ 



Proof. Assume r > p. Then 



1 /min(r, p) 
max(r, p) 



[cos k6 cos A;0±sin k9 sin /c0] . 
(A.5) 



ln(r^ + p^ - 2rp cos(6' =F </>)) = In r^ + In f 1 + - - 2- cos(6' =F (j))) 

\ r r / 

= 21nr + ln [(l - ^e^(^^'^)) (l - ^e^^^^^))" 

= 21nr-2^-('^) cos[A;(^T0)], (A.6) 

fc=i 

wherein the logarithms were expanded for p/r < 1. If instead r < p, we may simply 
interchange them in (A. 6). D 

Relation (2.18) is an immediate special case integral from Lemma 3, putting ^ = 0. 

More generally we have the following, putting Li2(-2) = Xl^i z'^/n'^ ~ ~ lo P^(-'^ ~ 
t)/t\dt (e.g., [2]). This dilogarithm function has a branch point at z = 1 and can be 
analytically continued through out the complex plane. 
Proposition 6. For r > p, 



/•TT 

/ ln(r^ + P^ - '^rp cos{9 =F (p))d(f) 
Jo 



27r\nr ±i 



Lh[-e-^'P)-Lh(e-^'P 



Lh[~e^'^P)+uJe'^'P 



(A.T) 



Proof. Using Lemma 3, we have 



ln(r2 + p2-2rpcos(^T0))rf0 = 27rlnr^2^^ 1^ ' ^ (^) sinfc^. (A.8) 



fc=i 



A;2 
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Then using the series definition of Li2 gives the Proposition. D 

In connection with (2.17), we have the integral 

U2{cp)pdp = — [4cVLi2(cp) + 2(cV - 1) ln(l - cp) - cp{2 + cp)]. {A.9) 

This follows from 

f^ " °° ^n ■n+2 

y 

^^ n^(n + 2) 

n=l ^ ' 

1 1 1 






lE'" 



2 

n=l 



p"+^ (A. 10) 



_n2 2n 2(n + 2)_ 
We also record the following result that is useful when r < p. 
Lemma 4. 

/ pLi2 la-\dp=-<. 2p^Li2 ( — j + ar[p + ar ln(p - ar)] - p^ In ( 1 j L 

(All) 
Proof. We use three integrations by parts: 

/ pLi2 ia-\ dp = -- / pln(l - ar / p)dp + — Li2 ( — j 

I [ arp P^ ^ I. I \ ^ P^j- f^^\ 

arp a^r^ n P^ , /-, / \ P^t • /a^A 

= — + ^ ln(p - ar) - — ln(l - ar/p) + y L12 ( — ) • 

D 
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